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^ ■ Abstract 

Q: 
o : 

^SJ , We consider the covariant Spectator equation with a One Boson Exchange 

' (OBE) kernel, and apply it to the NN system. Relativistic effects such as 

^ ' retardation and negative-energy state components are included in that equa- 

O . tion. We develop a method, based on the Pade method, to solve the Spectator 

(<~^ , equation without partial wave decomposition. The convergence of the partial 

I wave decomposition series is tested as a function of the energy. The on- and 

' off-mass-shell amplitudes are calculated. The NN interaction was fitted to 

' the differential cross section for NN scattering, in the energy range below the 

r~| . pion production threshold. 

^ ' 

^ I. INTRODUCTION 

c : 

^ ■ Electron scattering experiments off light nuclei at momentum transfer in the few GeV/c 

^ . region will provide new information on the intermediate and short range behavior of the 
nuclear interaction [1,2]. In this regime the interpretation of data requires a description 
based on relativistic kinematics and dynamics, and therefore sets an important challenge 
to theory. It also raises a technical computational difficulty, since traditional numerical 
methods to solve dynamical equations are based upon partial wave decomposition (PWD) 
of the NN interaction. In fact, scattering and bound states are often described by partial 
waves series. This procedure allows to explore nuclear symmetries and probes expansions 
in powers of energy /momentum. Moreover, by reducing the dimensionality of the integral 
equations it saves numerical computational effort, since for low energies (< 200 MeV) the 
series can be truncated to a few terms. However, for higher energies a large number of partial 
waves is required and the method becomes unfeasible, specially in applications for nuclear 
systems with A > 2 [3] . In fact two-body scattering amplitudes are input for the three-body 
amplitudes and, as it is shown in Ref. [4], above 300 MeV the convergence of the series 
requires more than 15 partial waves. Fortunately, present day computational limitations 
are not a serious objection anymore to obtaining directly the solution of three-dimensional 
integral equations without partial wave series. 
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In short, methods not based on PWD generate numerical problems which are i) straight- 
forward, presenting much less analytical and algebraic complexity than alternative formu- 
lations with PWD, ii) usable in the energy region where PWD converge too slowly, and 
which is probed by the Jlab data for electro- and photo- desintegration of the two- and 
three-nucleon systems, iii) solvable within reasonable CPU time in present day available 
computer power. 

Since for the energy range mentioned in point ii) relativity becomes also an issue, in this 
work we develop a method to solve relativistic quasi-potential equations, without PWD, for 
the NN system, in particular the Spectator equation [5,6] with a realistic NN interaction. 
The present paper prepares for later applications to photo- and electro- reactions at high 
energies. 

By construction, the Spectator quasi-potential equation incorporates relativistic effects 
such as retardation and negative-energy state components (p-spin states). The focus here 
is on obtaining a working method to solve it without PWD, proceeding directly through a 
three-dimensional integration. This is the first work that combines the two aspects, relativ- 
ity and three-dimensional numerical methods, in a application to the two-nuclcon system 
scattering problem. While we consider here the particular case of the Spectator equation, 
the method can be extended to other quasi-potential equations. The numerical complexity 
of the relativistic problem is twice as large as the one of the non-relativistic case. The re- 
sulting equation couples 8 channels corresponding to the different final helicities and p-spin 
cases. 

In Sec. II an overview of three-dimensional methods and of quasi-potential equations is 
presented. In Sec. Ill the detailed presentation of the Spectator equation for fcrmions is 
given. The NN interaction is also discussed, as well as the strong form factors required by 
the stabihty of the solution. In Sec. IV the method of solving the integral equation without 
partial wave decomposition is explained. In Sec. V results for NN on-mass-shell amplitudes, 
NN differential cross section and NN off-mass-shell amplitudes are shown and discussed. 
Finally, Sec. VII presents the conclusions. 

II. BACKGROUND 
A. Three-Dimensional (3D) Methods 

For energies not larger than 300 MeV, numerical methods without partial wave decompo- 
sition, also referred as Three Dimensional (3D) methods [4], were applied to non-relativistic 
NN scattering with spin and isospin variables, using the Bonn and Argonne interactions. 
The authors used 2-particle hehcity states and the role of the symmetries of the scattering 
amplitudes in that basis to advantageously reduce the size of the numerical 3D problem. 

Importantly, their results showed that in the vicinity of 300 MeV some of those NN 
amplitudes required already 16 partial waves to converge, making the partial wave method 
impracticable, specially for use in 3 nucleon calculations. Indeed, in applications to the 3 
nucleon system presented on Ref. [3] the PWD method was realized to be inappropriate 
for 3 nucleon processes at energies higher than 250 MeV. At those energies, both two- and 
three-nucleon scattering amplitudes were shown to exhibit a strong angular dependence in 
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the forward and backward scattering angles, preventing an efficient and reasonable descrip- 
tion in terms of a few partial waves. Other non-relativistic calculations without partial wave 
decomposition for the amplitudes of the two- and three- scalar particle systems were pub- 
hshed [7,8]. In Ref. [9] relativistic effects of several quasi-potential equations are compared 
using the 3D method for scalar particles. An advantage of the 3D calculation for scalar 
particles based on numeric linear system methods is that it works on a relatively large range 
of energies (300-800 MeV) with, almost, the same number of grid points [7,9]. 

B. Quasi-Potential (QP) equations 

Over the past years different relativistic formulations have been investigated for applica- 
tion to light nuclei, which can be classified in two major categories: Relativistic Hamiltonian 
Dynamics [10-12] and Relativistic Field Theories [5,13-17]. Wc will concentrate here on the 
second ones. They are based on covariant field theories with effective nucleon and meson 
degrees of freedom. The non-perturbative character of the nuclear interaction definitely 
requires an infinite sum of meson exchange diagrams, which can be summed effectively by 
means of integral equations, as the Bethe-Salpeter [13] 4-dimensional equation. 

In principle, the kernel of the Bethe-Salpeter equation should include all irreducible 
diagrams derived from a considered Lagrangian. In practice, models inevitably determine 
the truncation of the kernel expansion and only a restrict set of appropriate meson exchange 
diagrams is then kept. 

An usual approximation of the Bethe Salpeter equation consists on restricting the inte- 
gral kernel to a sum of One Boson (Meson) exchange diagrams, which is often called the 
Ladder Approximation. Nevertheless, the this approximation may be questioned. In fact, 
the one body limit is not recovered when one of the particle masses goes to infinity [5,18]. 
On the other hand, the crossed-box irreducible diagrams may be important, as confirmed 
by calculations within the Feynman-Schwinger formalism [19]. Although restricted only to 
the bound state of scalar particles, Ref. [20] showed that there are some 3-dimensional inte- 
gral equations that conveniently rearrange the series of ladder and crossed-ladder diagrams, 
yielding results equivalent to solving the BS equation beyond the ladder approximation. [19]. 
One of those 3-dimensional integral equations is the Spectator equation. 

The quasi-potential (QP) equations are integral equations that are obtained from the 
Bethe-Salpeter equation through a restriction on the energy variable yielding 3-dimensional 
integral equations. 

In particular, in the Spectator equation this restriction is motivated by important can- 
cellations between box and crossed-box amplitudes. These cancellations are such that their 
sum reduces to ladder diagrams of the mass pole of the heavier particle or, for equal interact- 
ing masses, to ladder diagrams of the mass poles of the two particles. For scalar interacting 
particles in the one-body limit the above cancellation is exact order by order [5,18]. When 
one of the masses is much larger than the other, the spectator equation in the ladder approx- 
imation alone gives practically the same result as the full Bethe-Salpeter equation, as shown 
in Ref. [21]. This outcome shows the efficiency gained in the rearrangement of the series. 
Furthermore, it has been shown that for nucleons (with equal masses) the non-vanishing 
difference between the sum of box and crossed-box diagrams and the spectator ladder di- 
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agrams can be effectively represented by one heavier boson exchanges [22]. The Spectator 
equation, with an appropriate OBE kernel, is therefore suitable for applications to the two- 
and three-nucleon systems, as successfully done in Refs. [2,6,23,24]. 



III. SPECTATOR EQUATION FOR FERMIONS 

In the original form of the Spectator equation [5] one of the two particles is considered on- 
mass-shcll in the intermediate state. To satisfy the Pauli Principle a new version consisting 
of a set of two coupled equations was presented in Ref. [6]. Since here we evaluate the 
scattering amplitude with both particles on-mass-shell in the initial state, the two equations 
reduce to a single one. 

Following Ref. [6] the scattering equation corresponding to the amplitude where particle 
1 is on-mass-shell in the initial and final states is written as 

Ma'P',ap{p',P; P) = Va'P',ap{p',P; P) (3.1) 

1 p k Tfh — 

-2 j T^^^oc'p',a^pAp\k]P)Goc^p^^oc2,li2{k]P)Ma2l32,oc(){k,P]P), 

ai,a2 ^ 



where 



Va'/3',a/?(p',P; P) = Va';3',a/?(P',P; P) + (-l)%a',a/?(-p', P; P), (3-2) 



is the anti-symmetrized interaction kernel. In the factor 6 = (—1)^ / is the total isospin of 
the 2-nucleon system. In the notation used m is the nucleon mass, p and p' are the relative 
initial and final 4-momenta, P the total 4-momentum. Also, the indices that represent the 
Dirac components of particles 1 and 2 are respectively a' and for the final state, and 
a and f3 for the initial state. E]^ is the on-mass-shell energy corresponding to the three 
momentum k: = vm^ + I?. We considered the cm. reference frame where P = {W, 0) 
{W is the total energy). When both particles are on-mass-shell the relative 3- momentum is 
represented by p and then W — 2Ep. 

In Eq. (3.1) Gaii3i,a2i32{k'i P) is the 2-nucleon Spectator propagator given by 

Ga,p,,a2p.Ak; P) = Aa,a2 {P/2 + k) G (P/2 - k) , (3.3) 

where ha^cfziki) is the particle 1 positive-energy projector and Gg^p^{k2) the particle 2 Dirac 
propagator. In Eq. (3.1) the asymmetry due to the propagator G is only apparent, since the 
anti-symmetrization is implemented by means of the kernel in Eq. (3.2). 

To write the scattering amplitude in terms of the helicity and p-spin we perform the 
change of basis similar to the one of the Ref. [6] : 

E <'a'(P'' A'l)<V(p', K)M^'P',ap{p\P\ PXaiP, Al)<^(p, A2), 
a, (3 
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where uf with j = 1, 2 are the asymptotic states of the particles (see Appendix A), Ai (A'^^) 
and A2 (A2) are the initial (final) helicities respectively for particles 1 and 2. The indices 
Pj and p'j with j = 1,2 express the initial and final p-spin states for particle j. A similar 
expression holds for V. 

The scattering amplitude in terms of the helicity and p-spin reads: 

Mtr^C^x. b'' P) = ^kIC^'L ip', P-^P) (3-5) 

- E j4^nx^.xSp'^k,p)9\k,^)^^^^^^ 



where 



Note that in the Spectator formalism one of the particles is always on-mass-shell in a positive- 
energy state and that is why in Eq. (3.5) there is always a positive p-spin state in the initial, 
intermediate and final states. The diagrammatic representation of the Eq. (3.5) is given in 
the Fig. 1. The crosses on the lines mean that the corresponding particles are on mass-shell 
with positive energy. 

The kernel V contains two terms, the direct term and the exchange term (see Eq. (3.2)), 
which are represented in Fig. 2. In this work we use an OBE kernel including vr, a, p 
and uo meson exchanges. For tt exchange we consider a mixture of pseudo-scalar (PS) and 
pseudo-vector (PV) coupUngs, defined as: 

K{p', k) = A,7' + (1 - A.)^4p^7', (3.8) 

where < A^ < 1 is the admixture parameter. 

The kernel parameters are presented in table I along with the cut-ofi^ parameters of the 
form factors (see Sec. IIIC). The ttIt^, rup and rrit^ parameters are fixed by the experimental 
values for physical mesons. The remaining 11 parameters were fitted to NN scattering 
data. We may anticipate at this point that the best results obtained favor the PV coupling 
(Att = 0). The model used for the kernel is very similar to the one of the Refs. [6,23]. Explicit 
expressions for V can be found in Appendix B. 

A. Prescription for the Exchange Kernel 

For the direct term the 4-momentum transfer q is given by 

= {p[ - hY (3.9) 

while for the exchange term by 

f = {p'i-k2f. (3.10) 
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The relation of these four-momenta with the relative four-momenta is given by 



P'i = (^p',p') (3.11) 

p'^ = (W-E^,,-p') (3.12) 

^i = (^k,k) (3.13) 

k2^{W-E^,-k). (3.14) 

When particle 1 is on-mass-shell with positive-energy the momentum transfer reads: 

-q^^{p'-kY-{E^,-Ey,f (3.15) 

-f = (p' + k)2 - {E^, + E^- Wf, (3.16) 



respectively for the direct and exchange terms of the kernel (see Fig. 2). Wc may then 
conclude that, while the direct term has no singularities in the meson propagators (we 
have always //^ — > 0), the exchange term may have the singularity corresponding to an 
on-mass-shell exchanged meson (//^ = (f ). The condition /x^ = indicates a real meson 
production from ofT-mass-shell nucleon states. However, since a real meson production 
process is not allowed below the pion production threshold, the singularity of V has no 
physical correspondence. It is simply a consequence of the anti-symmetrization of the kernel 
(see Eq.(3.2)). As shown by Gross et a/ this spurious singularity is cancelled by higher order 
diagrams [6]. 

In the numerical applications of the Spectator equations two prescriptions were consid- 
ered so far to deal with the spurious singularities [6]: 

• the "principal part prescription": the singularity is included but only the Principal 
Part of the integral is kept. This corresponds to the replacement 

I IJ? — q^ — ie ^ i jj? — qp i^iJ? — (p — ie ^^^^^ ^ 

• the "energy independent prescription": the momentum transfer is modified in order 
to remove the singularity. This amounts to the replacement 

that is, the momentum transfer is the same for the propagator in both direct and exchange 
diagrams, leaving the meson propagator denominator independent of the two-nucleon system 
energy W (see Eq. (3.15)). 

We now introduce the on-shell prescription, an alternative adopted by us in the present 
work. 

Eq. (3.16) can be rewritten as follows: 

= (p' + k)2 - {E^, - £;k)' -{W- 2E^,){W - 2E^). (3.17) 

The last term of this equation generates the spurious singularities mentioned above, and 
vanishes when both particles are on-mass-shell either in the initial state {W — 2Ek) or in 
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the final state {W — 2Ep/). We therefore define the "on-mass-shell prescription", by taking 
for —q^ the expression 

-f^{p' + \^f-{Ep,-Ey,f. (3.18) 

With this choice the exchange kernel V is consistent with the Feynman rules in the on-mass- 
shell limit. This is a physical argument in favor of the proposed prescription. Furthermore, 
the plus sign in Eq. (3.18) implies that when the direct term in a particular interaction 
is dominant in the forward direction, then the exchange term dominates in the backward 
direction, as it happens for interactions mediated by scalars. We point out that this property 
is not satisfied by the "energy independent prescription", where the exchange kernel q is 
replaced by the direct kernel momentum q. Note also that the "on-mass-shell prescription" 
is energy independent. 



B. Properties and Symmetries of the Amplitudes 

Using parity, time reversal and particle interchange symmetries we can decrease the 
number of independent amplitudes [25,26]. In particular, parity invariance reduces the 16 
helicity amplitudes with p' = +1, p = +1 states to 8 independent ones, according to: 



Mi= M++,++(p', u; p, 1) = M__,__(p', u; p, 1) (3.19) 

M2= M__,++(p', u; p, 1) = M++,__(p', u; p, 1) (3.20) 

M3= M+-,+-{p', u; p, 1) = M-+,-+{p', u; p, 1) (3.21) 

M4= M-+,+-{p', u; p, 1) = M+_,_+(p', u; p, 1) (3.22) 

M5= A^-+,++(p', u; p, 1) = (p', u; p, 1) (3.23) 

Me= M+-,++{p',u;p,l) = -M-+,^^{p'.u-pA) (3.24) 

Mr= M++,+-{p', u- p, 1) = -A1__,_+(p', u- p, 1) (3.25) 

M8= M-,+_(p', u- p, 1) = -A^++,_+(p', u- p, 1). (3.26) 

In the case p' = —1 the right hand side should include the phase (—1)^^ = — 1. 

The above relations are valid either p' is on- or off-mass-shell. Restricting to the on- 
mass-shell situation, further relations arise: 

Mj = Me (3.27) 

Ms = -Ms (3.28) 

Ms = -Mg. (3.29) 



The last identity is valid only for identical particles (exact isospin symmetry). As a result 
we are left with only 5 independent amplitudes. We point out that these relations are 
independent of the interaction used and of any prescription adopted, and are merely a result 
of the symmetries mentioned above. The relations (3.19)-(3.26) were tested numerically as 
a check to the code, since we did not explicitly impose the symmetries to reduce the number 
of equations. 
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If we consider instead a QP equation of the instantaneous type (no retardation), e.g. the 
Blanklcnbcclcr-Sugar [14] or the Equal-Time [15] equations, where the particles have always 
the same energy, we could use more symmetry properties to reduce further the number of 
off-mass-shell amplitudes. 

Finally we note that, for equations of the instantaneous type, a convenient combination 
of the helicity states defines states of well defined parity and two-body spin and helicity. The 
use of that basis states reduces the size of the numerical problem by block-diagonalizing the 
set of equations. This was done for example in Ref. [3,4] in a non-relativistic framework (no 
p' — —1 states). Nevertheless, for 3-body applications [3] the hehcity combination has to be 
inverted by calculating back the amplitudes in the asymptotic basis of uncoupled helicities. 



Nucleons are not elementary particles and their hadronic structure has to be taken into 
account through form factors. Mathematically, such form factors provide the necessary 
regularization of the integrals for the high order loops. Since in this work we solve the 
scattering equation without partial wave expansion, a careful study of the integrand function 
had to be performed in order to determine form factor functions adequate for the convergence 
of the method used. 

The starting point for the choice of form factors was the decomposition suggested by 
Riska and Gross [27] 



momentum, f^i is the form factor of meson i and /at the nucleon form factor. This factor- 
ization is adequate to describe the electromagnetic interaction with nucleons, and the scalar 
functions /at and could be interpreted as strong interaction vertex corrections and self- 
energy contributions of nucleons and mesons. As shown in Ref. [22] , however, some care has 
to be taken in the functional forms of the form factors considered in the above prescription. 
In fact, the most used momentum dependence of /at [6,23] induces additional spurious sin- 
gularities in a 4-dimensional framework, which can be eliminated by distorting conveniently 
the contour for energy integration around them [28]. These spurious singularities, however, 
do not occur within a QP context. 

For the meson form factor our choice is 



where Aj„. is the cut off of meson i. Note however, that by construction the functions frmio^) 
only modify the kernel for large (for 5^ = we have fmi{(f) — 1), and do not suppress 
the large momentum dependence of V{k, k; P). 
For the form factors we take 



C. Strong Form Factors 





(3.31) 



fN{k])= J 




(3.32) 



8 



with 



= - ml (3.33) 

being A^v the nucleon cut-off. The functions regulate the asymptotic behavior of both 
V{k, k; P) and V{k,p; P). The factorization (3.30) has been apphed by Gross and collabora- 
tors in several applications of the Spectator equation [6,23]. The function (3.32) with n = 1 
has been used for the first time in Ref. [23] and kept in Spectator equation applications since 
then. The meson form factor used in the same applications differs from (3.31) for large q^. 
The factorization (3.30) with the presented fm and /jv {n — 1) functional forms has been 
also used in Ref. [29], but in an instantaneous quasi-potential framework (where k = (0, k)). 
In the present calculation we found that in order to solve the Spectator equation with 3D 
methods the value n — 1, used in previous applications with PWD, was not large enough 
to allow the convergence of M.{k,p;P) for large |k|. We had to take instead n — 2. This 
is mostly due to the meson propagator behavior which peaks for forward and backward 
scattering angles at high values of the momentum [30]. Using the PWD method the meson 
propagator peak is smeared by the angle integration, but in the 3D method the structure of 
the propagator cannot be smoothened. 



IV. SOLUTION OF THE INTEGRAL EQUATION WITHOUT PARTIAL WAVE 

DECOMPOSITION 

In order to solve the scattering equation analytical or numerically we need to specify 
the scattering conditions, that is the initial and final momenta. We choose a reference 
frame where the incoming momentum is along the z axis and express the initial, final and 
intermediate momenta in terms of spherical coordinates 

p=(p,0,0) (4.1) 
p' = {p',9',^') (4.2) 
k=(k,^,^). (4.3) 

where p' = |p'|, p = |p| and k = |k|. 



A. Integration of the Azimuthal Angle 

To perform the (p integration we need to apply on the scattering amplitude a rotation of 
an angle ip around the z axis. In Appendix C we show that: 

Mt±XZi^, 9, p; W) = e^^'^-'^^'^Mt±XZi^, 9, 0; p; W), (4.4) 

where the p angles arc omitted by simplicity. 

Inserting Eq. (4.4) in the Spectator Eq. (3.5), and taking cp' = 0, we can factorize the cp 
dependence obtaining 
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Al^liS, (P'. 0; p; W) = V;*S,(p', 0; p; W) 



(4.5) 



p,A3,A4 ' 



k^dk f d |i^^M,H-.^^(p. k, 0; A, W^),^(k; VF)AiaS,,(k, ^, 0; p; W), 



where 



and 



y+P2'+P 
A^A2,A3A4 



(p', 9', k, e-\W)^^j dip e^~^^Vi$CZx,ip', k, ^, V'; (4.6) 

Ai — A2 



A 



(4.7) 



can take the values 0, ±1. 

The scattering equation, Eq. (4.5), is now in 2-dimensional form and includes the prop- 
agator function ^'''(k; W) which has a real pole at k = p [W = 2Ek) for p = +1. Performing 
the contour integration we obtain 

■M?I^A?A.(P', 0', 0; p; W) = V,?4';tX(p', 9', 0; p; W) 



■ E r 



AW 



A3,A4 



27r 

dcos 



2n 



V^yS'sX^iP'^ k, ^; A, W^)y.%(k; W^)Ma:,lU.(k, ^, 0; p; W) 



-V^v^ •tt,(p', P, e; A, W^)Mt3l;i?A.(P, 0; p; W). 



27r '^'lA2.-^3A4^ 



(4.8) 



In this equation we use g^^Q{k; W) to represent the £ = limit of g^'ik; W) (see Eqs. (3.6)- 
(3.7)). The multiplicative factor of the last term is a result of the residue 



1 / m \ ^ 

2 \E^) 



m^p 



(4.9) 



Jk=p 



In Eq. (4.8) only the function g^{k; W) has a singularity but we include the principal part 
symbol V also for p = — 1 for the sake of simplicity. 

The (^-integration (4.6) can be performed either analytical or numerically. In Refs. 
[3,4,7-9] the numerical integration was made. Here we choose the analytical integration in 
order to minimize computing time. The complexity of the analytical expressions depends 
crucially upon the fm form factors. Relatively simple results are obtained for the monopolar 
choice (3.31). In this case we can write Eq. (4.6) as 

d^e'^^Vl','f'l {p',p; P) = V Cr / dip- ^- r^, (4.10) 

^ AiA2,A3A4Vi^ 'i^' J r (a-6cos^)(c-6cos(/?)2' ^ ^ 

where F = 0, ±1, ±2, and Cp is a known function of the momentum magnitudes and polar 
angle. Any term of the last equation is subsequently easily integrated over. Meanwhile, in 
Ref. [31], where the equation for vr-N scattering is solved also without PWD, the authors 
came also to the same analytical technique but for an integrand not including form factors. 
Details of the analytical structure of the kernel can be found in Appendix B and the main 
steps for the analytical integration are given in Appendix D. 
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B. Numerical Method: Pade Method 



The relevant variables in the scattering equation (4.8) may be emphasized by the follow- 
ing convenient change of notation: we denote each of the eight possible helicity and p-spin 
combinations by a single index I' — {pg, A'^, A2}, Ik — {p, A3, A4} and Iq — {p2, Xi, A2}, ac- 
cording to table II (/', /fe and /q have nothing to do with the total isospin /); wc omit the 
total energy W, incoming momentum p and total isospin / dependences and perform the 
following substitutions 

VAt4'S.(p^^^o;p;w^)-v,,,„(p^«) 

y^!iXM^^'^ k, ^; A, W) ^ Vr,:,{^\ u; k, v), 

where 

u^cos0' (4.11) 

v^cosO. (4.12) 

Finally Eq. (4.8) simplifies to 

- E ^ j_^^Vi,j,{p',u;k,v)gH^;W)Mi,,i,{k,v) (4.13) 

m^p di 



-I- 



4^ 4^1 



-1 2^^^'A^P''^'P'^)-^^'^-fo^P'^)' 



with 



^ ^ ' ^ I ^ (k;VF) if 4 = 5,..,8 ^ ' 

In order to obtain a numerical solution of Eq. (4.13) we carry out a discretization of the 
integral variables, k G [0,H-cxd[ and v G [—1, 1], and use a gaussian quadrature integration 
technique. We choose a grid of A^p -|- 1 points for all momenta, p', k and p, and a grid 
of A^u + 1 points for the angular variables u and v. With this discretization procedure we 
transform the integral equation (4.8) into an algebraic set of equations: 

M = y + C-M, (4.15) 

where M and V are the matrix vectors M.i'jo{ki',Vji) and ViijQ{kii,Vji). The C matrix can 
be decomposed as C — A + B. For Ik — 1, .., 4, we have 

Ai'i'j'Wkij) = -j^kfVpj^ih, uy\ ki, Uj)g+{ki; W) (4.16) 
B{i'i'j'Uhij) = - ^27rV 2 VW ~ ^^''ikiki',Uj'; P, Uj)kNp+i (4-17) 



11 



with 



AS^S-S' (4.18) 

^ /■+°° „ ^ 1 1 , W -2m 
S^V dk^ = log (4.19) 

For Jfc = 5, .., 8, we have 

Ai'i'j'Wkij) = - J^kfVr,i,ih,Ujr, ki, Uj)g-{ki- W) (4.21) 
B{i'i'j'),{ikij) = 0- (4-22) 

In the previous equations w'^ and hj are respectively the gaussian weights for the variables 
ki and Uj. The momentum grid is obtained from a e ]0, 1[ grid by a change of variable 



i.- Xi 

where typically we take A ~ 0.5 m. In order to determine the on-mass-shell forward scat- 
tering amplitude we add the points kNp+i — P and mat^+i = 1 with zero weight. 

The dimension of the above matrices is n = 8(A^p + 1)(A^„ + 1), which is a large number 
when Np and Nu are of the order of 20. Therefore, a standard matrix inversion method 
requiring a large computer memory (for double precision complex numbers), becomes im- 
practicable. To overcome this limitation we use instead the Pade method, which gives a fast 
estimate of the result of the Born series for the coupled set of equations 

M = XV + XC ■ M, (4.23) 

where the A parameter is introduced by convenience and set to 1 at the end of the calculation. 
The usual power expansion for 2A^ -|- 1 terms reads 

M = AM(i) + A^M^^) + A^M(3) + ... + A^^+^M^^^+i). (4.24) 

The M^^^ vectors are evaluated by 

M(^) = V (4.25) 

M(^+i) = C-M(*). (4.26) 

and any element m of the M vector given by 

m = Ami + A^m2 + A^mg + ... + A^^+^msiv+i, (4.27) 

is approximated by a rational function of A 

/,x ao + Afli + ... + A^oat oqn 
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where the 2N + 1 a/ and hi coefficients are determined through the 2N + 1 coefficients, 
after equating Eqs. (4.27) and (4.28). This method is known in the hterature as SPA (Scalar 
Pade Approximant) [32] and mpade denotes the Padc[A^, iV] result. 

The advantage of the Pade method is that it replaces the matrix inversion by a fast 
estimate of the Born expansion, where all terms are evaluated as a matrix- vector multiplica- 
tion. This multiplication can be performed as n dot-products of two vectors of n dimension. 
Therefore the calculation requires memory to allocate only 2n, instead of n^, complex num- 
bers. This reduces substantially the dimension of the problem. The price to pay is the 
recalculation of the matrix lines. 

As we will see in the following section, 11 to 15 Pade terms are needed for convergence. 
The full calculation takes Ih to 2h CPU time in a Pentium IV at 3GHz. The memory size 
required is 50 MB. 

V. RESULTS AND DISCUSSION 

In this section we discuss the results obtained for the on- and ofT-mass-shell amplitudes, 
the fit to the np differential cross section, and the convergence of the amplitude as a function 
of the energy. 

The numerical results were checked to satisfy the Optical Theorem: 

MMtXZ^M 1; 1)] - E £ ^\Mt:>::,Ui^, v; p, i)p. (5.1) 

Since the optical theorem only probes the on-mass-shell amplitudes, the results for the off- 
mass-shell amplitudes were tested by numerically checking that they satisfy the symmetry 
properties as described in Section IIIB. 

A. On-mass-shell Amplitudes and Differential Cross Section 

Asymptotically the state of the two nucleons is characterized by the individual isospin 
states. Therefore, for the np system, we must consider 

nrA^,A,A.(p', 1) = In'y^Mx.ip' , p, 1) + lT'x'x',M,^ip', p, i), (5-2) 

where T^Py represents the anti-symmetrized matrix Ai with a total isospin I [Iz — 0). 

Figs. 3 and 4 show the scattering amplitudes (both real and imaginary parts) obtained, 
at a fixed energy Tiab = 300 MeV, for all independent helicity channels. The comparison 
between the exact result and the PWD results with increasing values of the total angular 
momentum J is shown. 

The convergence of the Pade amplitudes has been carefully examined, and the iteration 
procedure stopped when both real and imaginary parts of T converge with a relative error 
lower than 10~^. We conclude that in general grids with Np = 20 = 30 were enough for 
the numerical convergence of the solution, although Np = = 20 were sufficient below 200 
MeV. 
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For the helicity A1A2 = H — channel convergence requires 13 Fade terms for J = and 
11 terms for 1=1. For the A1A2 = ++ channel 15 terms for / = and 13 terms for / = 1 
are necessary. In all cases optical theorem has always satisfied with a relative error smaller 
than 10-2. 

After calculating all the on-mass-shell polarized amplitudes we can evaluate the differ- 
ential cross section 

In this equation IT"'?^ is given by 

Ai, A2 
Ai, A2 

which corresponds to an average over the initial helicity states and a sum of the final ones. In 
Fig. 10 the fit of the NN potential model to the differential cross section np data at energies 
Tiab—^^: 200 and 319 MeV is shown. Data are collected from the Nijmegen data-basis [33] 
and correspond to Refs. [34-36]. 

At least for the two first energies the quality of the fit is very good. The Tiab = 319 MeV 
energy case is already very near the pion production threshold, which justifies the slight 
decrease of the quality of the fit. Nevertheless, we conclude from this fit that the method 
developed in the present work is reliable and promising. Another point worth mentioning is 
that the fit selects the PV pion-nucleon coupling (mixing parameter Att = 0), in agreement 
with requirements from chiral symmetry. 

B. OfF-mass-shell Amplitudes 

The off-mass-shell T^^y AiA2(p'' ^' P' -*-) amplitudes may be plotted as functions of 2 vari- 
ables in momentum space. 

The results for the 8 p' = +1 amplitudes are presented in Figs. 6 and 7 for Tiab — 300 
MeV. The p' = — 1 amplitudes plots are presented in Figs. 8 and 9. Note that the on-mass- 
shell region corresponds to the line p' =0.375 GeV. 

It is interesting to note that the off-mass-shell region is very important for all channels. 
In fact the magnitudes of some amplitudes are important up to momenta of the order 1.5 
GeV, much larger than the on-mass-shell momentum. It is also important to point out that 
the amplitudes involving transitions to negative-energy states, p' = —1, have magnitudes of 
the same order of the p' = +1 amplitudes. Therefore, degrees of freedom corresponding to 
negative-energy states may have considerable weight within covariant effective field theories. 
Consequently, although explicitly absent in non relativistic models, they are accounted for in 
an effective way through the fitting procedures yielding low energy non relativistic realistic 
potentials. 
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VI. BREAKING OF PWD FOR HIGH ENERGIES 



We present here the study of the PWD convergence as a function of the energy. To 
perform this decomposition we follow Ref. [26]. Figure 10 shows the convergence of the 
PWD to the full calculation for a particular neutron-proton on-mass-shell T matrix, for 
three energy cases. We notice here that the imaginary part of T converges faster than the 
real part. For each energy case the criterium for convergence was defined as a deviation less 
than 1% from the full result. 

We can see in the Fig. 10 that for 200 MeV more that 10 partial waves are needed 
to achieve convergence. We confirmed that some hehcity channels at 300 MeV require 16 
partial waves, as obtained in the non-relativistic calculation of Ref. [4]. As mentioned before, 
this large number of partial waves makes the PWD method impracticable to generate the 
two-body amplitudes required as an input for the three-body bound or scattering calculation 
as seen in Ref. [3] within a non relativistic framework. 

On the other hand, the results for 100 MeV show that it is a good approximation at low 
energies to take J = 6 as the highest total angular momentum. This explicit finding fully 
justifies the cut-off at J = 6 reported in the three-body bound state relativistic calculations 
of Ref. [23]. 

VII. CONCLUSIONS 

We considered the Spectator equation for nucleons. Its solution depends on the helicities 
of the particles, as well as on their p-spins. {p' — for the positive- energy state component 
and p' = —1 for the negative- energy state component). 

The scattering amplitude is solved without partial wave decomposition using the Pade 
method. This method revealed to be efficient and suitable for the solution of the relativistic 
integral equation without partial wave decomposition. For Tiab < 350 MeV, from 11 to 15 
Pade terms were needed for convergence. 

Strong form factors and a prescription for the exchange kernel different from the ones 
considered in other calculations within the spectator formalism are used. When both parti- 
cles are on-mass-shell the present prescription for the exchange kernel gives rise to the kernel 
directly obtained from the Feynman rules. This is important in view of possible applications 
to problems where the fields are not effective, the qq bound state and the quark exchange 
diagrams in vrvr scattering. 

We fitted an OBE interaction to the differential np differential cross section in the 100- 
350 MeV energy range. We let the percentage of PS and PV admixture to float as a free 
parameter of the fit. It turned out that the PV coupling was favored by the fit, which is in 
agreement arguments of chiral symmetry. 

We studied the convergence of the PWD method as a function of the energy. We con- 
cluded that above Tiab = 250 MeV at least 16 partial waves have to be included for some 
helicity cases. This large number of terms indicates the breaking of the PWD method for 
applications in heavier nuclei. 

Beyond the p' — +1 amplitudes involved in the cross section calculation, we have also 
calculated the p' — —1 amphtudes related to processes involving one off-mass-shell particle. 
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These amplitudes are numerically significant and may be used in the calculation of meson 
production cross sections [37,38] and of the ^He(e, e')X reaction observables for 1-10 GeV 
energies [1]. 
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APPENDIX A: STATES OF HELICITIES 
1. Positive-energy state spinors 

Following the construction of Refs. [39,40] we obtain for the spinors the expressions: 



Uj(k, A) = Afe 



1 



|A >j, (Al) 



with the normalization 



The Pauh spinors of the particle 1 and 2 

|A >i= XA(k) (A4) 
|A>2=^A(k). (A5) 

are related by 

V'A(k)=X-A(k). (A6) 
Initial and final state Pauli spinors are presented in table III. 



2. Negative-energy spinor states 



The negative-energy spinors are constructed from positive-energy states using the charge 
conjugation operator C [6,39,40]: 

^;i(k,A) = (-l)AC4(k,A) (A7) 
V2(k,A) = AC«f(k,A), (AS) 
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where C — — 27^7^ and T indicates matrix transposition. The relative factors are introduced 
by convenience. As a result we get 



(k, A) = iVfe 



1 



lA >i 



(A9) 



APPENDIX B: KERNEL V 



In this Appendix we describe the analytical structure of the kernel with the OBE form. 
As mentioned before, for any isospin state /, the kernel V contains the direct, V, term and 
the exchange kernel, V, term: 



For meson i the direct term is given by 



;+ p',+ P 

where the exchange momentum is 

= {E^, - E^f - (p' - k)^ 
For vector mesons we have to undertake the substitution 

{p'l - - k2)u 



A1A2 ^ Klip'^, ki)A'M, k2) 9i.u + 
Similarly for the exchange term we have 



u^ip'i, \[)Ai{p[, k2)u^{k2, A2) 



U^'{p'„ X',)A2{p'2, k,)u+{k,, Xl)[fmAfWfN{p'2')fN{kj), 

where, according to the " on- mass-shell prescription" the exchange momentum is 

= {E^, - E^f - (p' + k)2. 
For vector mesons the replacement is now 

(pi - k2)^{p'2 - ki)u 



A,A2^A'i{p\M)A''2{p'2M) 



(Bl) 



(B2) 



(B3) 



(B4) 



(B5) 



(B6) 



(B7) 



The calculation results are presented in Sec. B 1 for the direct term and Sec. B 2 for the 
exchange term. For the sake of simplicity we do not explicit the parameters involved but 
only the analytical structure. 
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1. Direct Kernel 



We will write the direct kernel by means of the following auxiliary functions 

^?(p',k) = x'4(p')XA,(k) (B8) 
Z°(p',k)=V'y;(pOV'A.(k), (B9) 



^i(p',k)=x'4(p')<Tf^XA,(k) (BIO) 

Z^(p',k)=^4(p')afVA.(k) (Bll) 

The functions (j = 1,2, a = 0, ..,3) are calculated from the Pauli spinors presented in 
Appendix A. For simplicity sometimes we suppress the arguments of Z^. 

a. Meson a 

Vt'SiMMip'^k) = --^,[U{q')ffN{j>'')fN{e) 
l-''a y 

Nl,NlH„{p\k) (B12) 
Z°(p',k)ZO(p',k) 



Hfj{p', k) is also a function of A'^^, A2, Ai, A2, p' and p. 
For large k we have H„{p' , /c) ~ 1. 



h. Meson tt 
•,2 



N}N'Mp',k) (B13) 
Z?(p',k)Z°(p',k) 



Ht^{p', k) is also a function of A'^^, A2, Ai, A2, p' and p. 

For large k we have H„{p', k) ^ 1 for PS coupling and i?7r(p', k) ~ k/m for PV coupling. 



c. Vector mesons 

Nl,NlH,{p\k). (B14) 

i=l j=l 

+r3 E(K + ^i)^i^2° + r, jZ{p'i + ^i)'^?^2°- (B15) 



i=l 
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The coefficients r/ {I — 0, ..,4) are functions of the momenta p', k, of the indices A'^, Ag, 
'^i) -^2) p') P and of the parameter. 

For large k we have Hy{p', k) ^ 1 when Ky = and Hy{k, k) ~ k^/m^ and Hy{p', k) ~ k/m 
when Ky ^ 

2. Exchange Kernel 

For the exchange term we use the auxihary functions 

Z°(p',k) = x'4(p')VA.(k) (B16) 
Z°(p',k) = V^;|(p')XA,(k), (B17) 

Zi(p',k) = x'4(p')^fVA.(k) (B18) 
Z^(p',k) = ^4(p')af)xA.(k). (B19) 

These functions are also evaluated from the Pauli spinors. 



a. Meson a 



N^,N^H,{p',k) (B20) 
ZO(p',k)ZO(p',k) 



H„{p', k) is also a function of A'^^, Ag, Ai, A2, p' and p. 
For large k we have H„{p\ /c) ~ 1. 



h. Meson vr 

P-K y 

N}Ni%{p',k) (B21) 
i°(p',k)i°(p',k) 

Hn{p', k) is also a function of A'^, A2, Ai, A2, p' and p. 

For large k we have H.^{p' , A;) ~ 1 for PS coupling and H^^ip', k) ~ k/m and H^^lk, k) ~ 
k^/m^ for PV coupling. 
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c. Vector mesons 



9v 



N}NlH,ip',k). 



(B22) 



i=l i=l 

+r3 E(p: - fc^)^i^2° + ^4 jZip'r - hfZlZl (B23) 

i=l i=l 

The coefficients r; (/ = 0, ..,4) are a function tiie momenta p', k, of A'j^, Ag, Ai, A2, p', 
p and also of the parameter. The r; coefficients for the exchange kernel are not related 
with the direct kernel coefficients. 

For large k we have Hy{p' , /c) ~ 1 when Ky — Q and H„{k, k) ~ k^/m^ and Hy{p' , k) ^ k/m 
when Ky ^ 0. 



APPENDIX C: (p ROTATED AMPLITUDE 

In this appendix we derive the relation between the scattering amplitude on the scattering 
plane (p is on the z axis) 

and the scattering amplitude on a rotated plane characterized by (/^'-rotation in the z axis 

Thus, we consider the Lorentz transformation 

A = (CI) 

corresponding to a rotation of an angle ip' around the z axis. 

The correspondence between the spinors before and after this Lorentz transformation is 

S{A)u1{p, A) = E D,,,(i?A)<(Ap, A') (C2) 

A' 

^(A)«^(p, A) = J2D-x',-x{RaH{Ap, A'), (C3) 

A' 

where S{A) is the operator that transforms the u, v-states in the A Lorentz transformation 
and D is the usual D^/^ Wigner matrix in terms of the rotation angles. The rotation 
operators are given by 

i?A = Hj[^AHp (C4) 
R'^ = Hxi,AH,,. (C5) 
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In the last equations Hp is the operator that transforms a 4-momentum (m, 0) into p — 
(i?p,p). Details can be found in Ref. [41]. The operation can always be written in a 
sequence of a boost (L^) and a rotation (i?p) 

Hp — RpLp. (C6) 

After the Lorentz transformation has been done the following relation between the orig- 
inal and the rotated scattering amplitudes is obtained 

Ma'p'M^p',Ap;AP)= (C7) 

Oil, Pi 
0:2, P2 

for more details see and Appendix B of Ref. [6]. 

Changing (C7) to the helicity representation according to eq. (3.4) and using the operator 
rotation properties and the invariance of permutation between boost and rotation operators 
related to the same axis, we finally conclude that 

Mt/^Sliv', 9', p; W) = e^^-'A^+^t - (p', 9', 0; p; W). (C8) 



APPENDIX D: FUNCTION V 



In this Appendix we describe how to evaluate V defined by Eq. (4.6). By performing 
the (fi integration we get: 

V{p', 6', k, 6; X, W) = r e^~'^V{ip)dip, (Dl) 

^0 

where we use the simplification 

V(^) = VaT^'JI^p'' 0; d, f, W). (D2) 

We note that 

X = (D3) 

SO A = 0,±1. 

First, we separate the y9-dependent parts from the independent ones (factorization). 
Next, we analyze the structure of the resulting functions. We need to consider two different 
cases: the non vector meson exchange {a and tt) and the vector meson exchange. 
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1. Factorization of V{(p) 
From the V expression (see Appendix B) we conclude that 

V(^) = 5iAtgl,V(^)n, (D4) 

where 

m - ^^""if^^^ (D5) 

and 

n = N}N',Mp")Uk')H,ip', k). (D6) 
In this equations i is the meson index. Also from Appendix B we have 

Hi{p',k) = H,{p',k)ZlZl (D7) 

for a non vector meson and 

Hi{p',k)=H,{p',k), (D8) 

for vector mesons. Note that the parameterization (D6) is valid for the direct term and for 
the exchange term if we replace by (f, Hi by (and by Z° for the non vector case). 

Using the coordinates definition (4.1)-(4.3) with (p' — for and we can conclude 
that 

V{v) = ^ 7 ^, (D9) 

a — b cos (p{c — b cos (py 

where 

a = ^ p' 2 + ip 2p'k cos e' cos ^ - (DIO) 
6 = ±2p'ksine'sine. (Dll) 
c = A^ + p' ^ + k^ If 2p'k cos e' cos ^ - g^. (D12) 

The upper sign should be used in the direct term and the lower sign in the exchange term. 

2. 7^ in terms of ip 

The next step is to write TZ of eq. (D6) in terms of (p. We need to consider two separate 
cases: the non vector mesons and the vector mesons. 
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a. Non vector mesons 

For non vector mesons we can write 

7^ = /(p^k)•Z?Z°. (D13) 

The exact expression of /(p',k) can be easily deduced from (D6). Attending to the Z^Z2 
dependence in (p, we conclude that 

n = /(p', k) • (Co + ci e-*^^^ + C2 e'^'^), (D14) 

where q (/ = 0, ..,2) are known coefficients depending on the scattering conditions and on 
the hehcities states. We can write 

V(p',e',k,e;~X,W)^5jAiglfip',k). 

[coJ^o(A) + ciJ^o(A - A3) + C2J^o(A + A3)], (D15) 

where the function J-o{n) is defined as 

1-27T 

J^o(n) = / ci<^e^"'^ V{^), (D16) 
Jo 

and n are n — 0, ±1, ±2. 

b. Vector mesons 

For vector mesons the function TZ can be written as a linear combination of the terms 

Z^'Z^^ and hZ^'Z^^ 

where ai,a2 = 0,..,3. The first term is reduced to the non vector case discussed in subsection 
D2a, because Z'^^Z^'^ can also be written as 

Co + Ci 6-*^=^^ + C2 e*^^'^, 

with appropriated coefficients. The second term can be decomposed in 3 cases considered 
as follows: 

Case 1 (/ci = k sin d cos 99) 

In this case we need to integrate factors hke 

(ksin^) cos99e^"^l7(<^), 

and the corresponding term of V ^ which we label Vi, is 

(p', e', k, ^; A, ) = bikUl.f{^, k) (k sin Q) ■ 

[co.Fi(A) + ci^i(A - A3) + C2^i(A + A3)], (D17) 
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where 

Tiin) = r dip cos LP e'^^^Viip), (D18) 

with n = 0, ±1 ± 2. 
Case 2 (A;2 = k sin 9 sin Lp) 

We have also to consider terms hke 

(ksin^)sin(/je^"'^y((^), 
from which results for the corresponding V, which we label V2 
V2{v\ e', k, 6; A, W) = 6jAiglf{p', k) (k sin 6) ■ 

[CoJ^2(A) + CiJ^2(A - A3) + C2J^2(A + A3)]. (D19) 

The J^2 function is defined as 

J^2(n)^ d(p sin (fie''''^V((p), (D20) 
Jo 

with n = 0,±1,±2. 
Case 3 {k^ — cos 9) 

No new (^-dependence appears. This case reduces to the non vector meson case. 

3. Functions Ti{n) 
The functions Ti{n) with / = 0, 1, 2 can be written in terms of the integrals 

Ri — '^i^) cos' (fdif, 
Jo 



for / = 0, ..,3. The Ri integrals are performed analytically with the software program 
Mathematica and simplified afterwards. 
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TABLES 





138 MeV 




13.470 


\ 


0.0 




1190 MeV 
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TABLE I. Model parameters. 
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TABLE II. Table defining the indices /q, h or 
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Initial state 

A = l A = -l 



Xa(z) 

-0^(2) 



I) 



Final state 

A' = 1 A' = -1 



cos |e 2 \ / — sin|e 2 

sin |e*2 / \ cos |e*2 



— sin '2 \ / cos 



V'A'l^,^) .^Jj^ ..Jej^ 



cos |e 2 y y sm |e 2 



TABLE III. Pauli helicity states. The phase differences between these spinors and the ones of 
these Refs. [39,40] are due to the fact that we consider, for convenience, the rotations convention 
of Ref . [40] and the phase convention for particle 2 of Ref. [39] . 
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FIGURES 

FIG. 1. Helicity representation of the Spectator equation. 
FIG. 2. a) Direct term of the kernel V. b) Exchange term of the kernel V. 
FIG. 3. Helicity amplitudes for 300 MeV and partial wave decomposition. 
FIG. 4. Helicity amplitudes for 300 MeV and partial wave decomposition. 
FIG. 5. Differential cross section results for 99, 200 and 319 MeV. 
FIG. 6. Off-mass-shell amplitudes for np process with p = +1, p' = +1. 
FIG. 7. Off-mass-shell amplitudes for np process with p = +1, p' = +1. 
FIG. 8. Off-mass-shell amplitudes for np process with p = +1, p' = —1. 
FIG. 9. Off-mass-shell amplitudes for np process with p = +1, p' = —1. 
FIG. 10. Partial wave decomposition of the M3 amplitude for 100, 200 and 300 MeV. 
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FIG. 1. Helicity representation of the Spectator equation. 
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FIG. 2. a) Direct term of the kernel V. b) Exchange term of the kernel V. 



1 



Full calculation 

J=4 

J=6 




T(++ -> ++) 



-0.5 

COS 



Full calculation 

J=6 

J=10 

J=14 



0.5 



T(++ -> ~) 



-0.5 0.5 

COS 9 



20 



> 



E -20 



-40 



-1 



10 



> 

(D 



-5 'V 



* Full calculation 
J=4 
J=6 



V. 



T(++ -> ++) 



-0.5 





COS 9 



Full calculation 



0.5 



J=6 
J=10 



/ 



/ 



/ 



/ 



T(++ -> ~) 



-10 



-1 -0.5 0.5 

COS 9 




-0.5 0.5 

COS 9 

FIG. 3. Helicity amplitudes for 300 MeV and partial wave decomposition. 
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This figure "fig6.gif" is available in "gif" format. 
FIG. 6. Off-mass-shell amplitudes for np process with p = +1, p' = +1. 



This figure "fig7.gif" is available in "gif" format. 
FIG. 7. Off-mass-shell amplitudes for np process with p = +1, p' = +1. 
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FIG. 8. OfF-mass-shell amplitudes for np process with p = +1, p' = —1. 

This figure "fig9.gif" is available in "gif" format. 
FIG. 9. Off-mass-shell amplitudes for np process with p = -|-l, p' = — 1. 
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FIG. 10. Partial wave decomposition of the M3 amplitude for 100, 200 and 300 MeV. 
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